Localization in discontinuous quantum systems 
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Classical and quantum properties of a discontinuous perturbed twist map are investigated. Dif- 
ferent classical diffusive regimes, quasilinear and slow respectively, are observed. The regime of 
slow classical diffusion gives rise to two distinct quantal regimes, one marked by dynamical localiza- 
tion, the other by quasi-integrable localization due to classical Cantori. In both cases the resulting 
quantum stationary distributions are algebraically localized. 

PACS numbers: 05.45. +b 



A major feature of quantum dynamics of classically 
chaotic systems is the quantum suppression of classical 
chaotic excitations, a phenomenon known as dynamical 
localization. A prototype model, both for classical chaos 
and quantum dynamical localization is the Kicked Ro- 
tator Model E| (KRM) , whose dynamics is described 
by the well known Chirikov standard map || (CSM). 
This is a 2-d continuous perturbed twist map, with a 
transition point, discriminating between bounded mo- 
tion (prevalently regular on invariant KAM tori) and un- 
bounded and diffusive one (prevalently chaotic). Even 
though transport properties of 2-d maps are now quite 
well understood, analytical results are only possible for 
particular maps, e.g. linear ||. In particular, the latter 
are the simplest discontinuous perturbed twist maps on 
the cylinder. For such discontinuous maps the hypoth- 
esis of KAM theorem are not satisfied and the motion 
is typically unbounded even if it is possible to mark two 
different dynamical regimes (both diffusive). Discontin- 
uous maps also emerge from the study of more concrete 
physical models, such as the motion of a particle collid- 
ing elastically within a two-dimensional bounded region 
(billiard Q). On the other side very little is known about 
the quantum dynamics of such discontinuous maps. In 
particular, it is far from being obvious that the relation 
between quantum localization and classical diffusion, ob- 
tained for the KRM, holds in this case too. 

To answer the above questions, let us consider the 
following discontinuous map on the cylinder [0, 27t) x 
[-co, oo] 



p = p + kf(0) 

9 = 9 + Tp mod - 2tt 



(1) 



where f(9) = sm(6) sgn(cos0). This function is a 
particularly simple approximation of the stadium map 
||,[|. Moreover, it is quite similar to the CSM (where 
f(9) — sin#) which has been widely investigated in the 
past. 

Even if the following analysis has been put forward for 



this specific function, it can be generalized @] to generic 
discontinuous, periodic and bounded < 1) func- 

tions. This set of functions can be also enlarged to con- 
tinuous bounded functions with a discontinuous deriva- 
tive. In this case the situation is slightly complicated, 
since usually a critical value of the parameter K = kT 
appears (see for the piecewise linear map) such that, 
when K < K cr , the phase space is covered by invariant 
tori which do not permit unbounded motion along the 
cylinder: only for K > K cr the motion is diffusive. 
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FIG. 1. Poincare surface of section for k — 0.01, T = 1. 
a) 10 different particles with initial momentum p — 0.001 
and different phase 9 have been iterated n = 10 3 times for 
the Chirikov standard map. (b) One particle starting at the 
point po — 0.011 and 6q = 3 have been iterated 3 • 10 4 times 
using the discontinuous function f(8). 

The CSM is characterized by unbounded diffusive mo- 
tion in the momentum p for kT > 1 while, when kT < 1, 
the motion is prevalently regular with regions of stochas- 
ticity bounded by KAM invariant circles. On the other 
hand, the classical properties of map (Q) are quite dif- 
ferent. Indeed, due to the discontinuities of f(6) at 
9 = vt/2,3/2?t, the hypothesis of KAM theorem are not 
satisfied and, generally speaking, KAM tori does not ex- 
ist, even for very small k. This means that one trajectory 
fills, in a dense way, any portion of the cylinder (phase 
space), for any k ^ 0. Nevertheless the phase space is 
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filled by Cantori (remnants of KAM invariant tori) 
which constitute partial barriers to the motion Due 
to the sticking of trajectories along these invariant struc- 
tures the diffusive motion is slowed down in close analogy 
to the saw-tooth map case described in [0 . 

An example of the classical map dynamics is given in 
Fig.l. In the right picture (b), the Poincare surface of 
section is shown for the discontinuous map (0). A sin- 
gle initial condition has been iterated n — 3 • 10 4 times. 
As one can see a single particle is free to wander in the 
whole phase space but the motion is far from being ran- 
dom. Indeed, due to sticking in the neighborhood of 
Cantori, the trajectory is almost regular on a finite time 
scale t. Diffusive motion results from jumping among 
different stable varieties belonging to different Cantori. 
As the iteration time, or the number of initial particles, 
is increased, regular structures disappear and the surface 
of section appears to be covered uniformly. For sake of 
comparison in the left picture (a) the same portion of 
phase space is shown for the CSM, with the same value 
of fc. Here 10 different trajectories have been iterated 
n = 10 3 times: each trajectory covers just one torus. 

Despite the "quasi" regularity of the motion, numerical 
results show that, when kT < 1, the dynamics is diffu- 
sive, for t > t, along the cylinder axis (p— coordinate) 
with a diffusion rate D given by 



D = li 



(P 2 (n)} 



= D k 5/2 VT 



(2) 



where n is the time measured in iterations of the map 
(|l|) and the average (...) has been performed over an ini- 
tial ensemble of particles with the same momentum p and 
random phases < 6 < it. Also, in (]|), Dq ~ 0.4 is a nu- 
merical constant (dependent from the function f(0)) and 
the factor y/T has been added for dimensional reasons. 

On the other side, when kT > 1, the random phase 
approximation Q can be applied and one finds diffu- 
sive motion along the p— direction with a diffusion rate 
D ~ D q i = fc 2 /2, where D q i is the diffusion rate in the 
quasi-linear approximation, namely assuming the phases 
8 to be completely random uncorrelated variables. No- 
tice that, in the undercritical region kT < 1, the diffusion 
coefficient D ~ k 2 \JkT is less than the quasilinear one 
Dqi ~ fc 2 , due to the sticking of trajectories close to Can- 
tori. In Fig. 2 the dependence of the diffusion rate D is 
shown as a function of fc for T = 1. The dashed and 
full lines indicate respectively the quasilinear diffusion 
(kT > 1) and the slow diffusion (kT < 1). 

The apparently strange dependence of D on fc, in the 
"slow" diffusive case kT < 1 was found in similar discon- 
tinuous maps, e.g. the saw-tooth map || (f(0) = 0/2n), 
or the Stadium map Q. In Ref. || a theoretical ex- 
planation of the exponent 5/2 was given in terms of a 
Markovian model of transport based on the partition of 
phase space into resonances. 
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FIG. 2. Diffusion rate for the discontinuous map as a func- 
tion of k and T = 1. Open and full circles indicate respectively 
the "quasilinear" and the "slow" diffusion. Dashed line repre- 
sents the quasilinear approximation D = D q i — 0.5fc 2 which 
holds for k> 1. Full line is the best fit D = 0.4fc 5/2 obtained 
from full circles. 

Let us now consider the quantized version of map (Q). 
According to a well-known procedure Jlj the quantum 
dynamics can be studied by iterating the quantum evolu- 
tion operator over one period Wt, starting from an initial 
state ipo(O) 



iP(T) = U T ^o 



-ihTn 2 /2 -ikV(0)/h 



to 



(3) 



In (||), as usual, n = —ihd/dO and V(6) = |cos(9|. 
Quantum dynamics depends on both parameters k/fi and 
Th separately. These parameters can be renormalized by 
letting k/h — > fc and Th — > T (which is the same as to 
put h = 1). The semiclassical limit is then recovered by 
performing simultaneously the limits fc — > oo and T — > 
keeping kT = const. 

The most studied example of quantization of twist 
maps like (0) is the KRM §, where V(6) = cos 8. Nev- 
ertheless the regime kT < 1 , differently from the case 
kT > 1, fc 3> 1, was not object of intense investigations. 

At least numerically, one can observe two different 
regimes, distinguished by the so-called Shuryak border 
fc = T H . For fc > T the quantum steady state is expo- 
nentially localized over a number l a ~ y/k/T of momen- 
tum states [^O|,0. This number has been interpreted 
jnH , in a realistic way, as the number of quantized mo- 
mentum states contained in the main classical resonance 
(see Fig.la) the size of which is y/k/T 0. For fc < T the 
width of the principal resonance is smaller than the dis- 
tance among quantized momentum levels, and no kind of 



2 



semiclassical excitation process, based on the overlapping 
of resonances, is possible. In the following, the analysis 
will then be restricted to the case k > T only. 

Before studying the discontinuous case, let us recall 
a few important facts related to the evolution operator 
. Due to the discontinuity in the first derivative of the 
potential V(9), the matrix elements oIUt in the momen- 
tum basis decay according to a power law away from the 
principal diagonal: |W„, n+s | ~ 1/s 2 . This case was inves- 
tigated |L2] for Band Random matrices: it was found to 
be typically characterized by power-law localized eigen- 
states around their centers no, |0 n | — I" — "o|~ 2 - 

The following question is then important : is it possi- 
ble to connect quantum localization lengths and classical 
diffusion rates, as in the case of the KRM? If so, what 
is the critical border necessary to start the classical-like 
diffusion process? Moreover, what is the role played by 
classical invariant structures, such as Cantori, in quan- 
tum dynamics? 
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FIG. 3. Localization length la as a function of y/k/T for 
fixed T — 0.0 1. Lines are the theoretical predictions : dotted 
(la = y/k/T), full (l a = D). Dashed line is the quasilinear 
border k q {T = 1. 

To answer the last question, let me remind the pioneer- 
ing works p3|jr^ ] where quantum propagation of wave 
packets through the classical Cantori was first investi- 
gated. Other important results can be found in JL5[ where 
it was proposed that Cantori could act, in Quantum Me- 
chanics, as total barriers to the motion if the flux ex- 
changed through turnstiles is less than fi. One can then 
reasonably assume that, in the deep quantum regime, the 
system will not be able to "see" the holes in the Cantori 
which behave as classical invariant tori. 

A more refined analysis requires the introduction of 
some kind of measure of the quantum distribution width. 
Since in this model localization is presumably not ex- 
ponential, a unique scale of localization is not properly 



defined. For instance, while in the case of exponential 
localization the usual measures of localization, e.g. in- 
verse participation ratio, variance, entropy pH , coincide, 
for power law localized distributions the dependence on 
the parameters can be different if different definitions are 
adopted. Then we choose the variance as a measure of 
the distribution extension (degree of localization) : 

t=E» a i^(*)r] 1/a w 

n 

which has a proper semiclassical limit. Since this is, in 
general, an oscillatory function of the iteration time, a 
further average in time is necessary in order to get time- 
independent results. 

Numerical data are presented in Fig. 3 where l a has 
been plotted as a function of y/k/T. Excluding oscilla- 
tions, data follow, for k < k cr , the dotted line y/k/T, 
as for KRM. Indeed, as one can see comparing Fig. la 
and Fig. lb the principal resonance and the "quasi" prin- 
cipal resonance have roughly the same size. This is a 
manifestation of the regularity imposed by quantum me- 
chanics, or, in other words, of the discrete nature of the 
quantum phase space. This means that classical discon- 
tinuous structures behave exactly as continuous ones. 

On the other side, since the classical discontinuous sys- 
tem is diffusive, the number of occupied quantum states 
should increase on going into the semiclassical region. 
Following known arguments for the dynamical localiza- 
tion, one can expect the localization length to be given by 
the number of states inside a "quasi" principal resonance 
(y/k/T), as soon as it equals numerically the classical 
diffusion coefficient. In this way the critical value k cr 
can be obtained by equating the following expressions: 

l a ~ y/k/T ~D = D k^ 2 VT (5) 

which gives the value k cr = \/ y/D$T . 

It is important to notice that the "quasi-integrable" 
value I u ~ y/k/T can survive well above the thresh- 
old k — 1 which is the value necessary to start the 
classical-like diffusion process for the KRM when kT > 
1. Also, this kind of localization is not connected with 
any classical-like diffusive process, resulting instead from 
a quasi-periodic motion. The absence of diffusive quan- 
tum motion, in the region T < k < k cr can be ascribed 
to a "dynamical " diffusion rate l a less than the size 
of the "quasi" -principal resonance l a ~ y/k/T. For 
instance, numerical simulation indicates a localization 
length l a ~ 80±10 ~ y/k/T for k = 10 > 1, T = 1/1000, 
while D ~ 0.4fc 5 / 2 v / T = 4. In other words, in the region 
dominated by slow diffusion, the threshold for classical- 
like diffusion is k > k cr — 1/ 'y/D^T and not k > 1. 

These theoretical predictions are confirmed by the nu- 
merical data presented in Fig. 3, which closely follow the 
curve (full line) l a = D for k cr < k < k q i- Here k q \ 
stands for the border of validity of quasilinear diffusion 
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: k q i = 1/T. This confirms and extends the validity of 
the dynamical localization theory even in the presence 
of "slow" diffusion and algebraic decay. This last point 
can be directly observed in Fig. 4a where the quantum 
steady state distribution P(n), is shown together with 
the corresponding line n~ 4 . 
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FIG. 4. Quantum stationary distribution for k = 50, and 
T = 0.01. The map has been iterated 10 6 times. The final 
distribution is obtained by averaging over the last 10 5 kicks. 
The initial state is ip n — 8 n ,o- (a) Log- log scale. The line n -4 
has been drawn to guide the eye. (b) Log scale. 



At the critical point fc cr , separating these regimes, quan- 
tum dynamics starts to follow the classical excitation pro- 
cess. Differently from the KRM, for which k cr ~ 1, one 
finds here k cr ~ 1/y/T. 

During the completion of this manuscript I became 
aware of another related work |l6| where a regime of 
quantum integrability is found, for the Stadium billiard, 
in the region delimited by the inequalities Ee > 1 and 
\f~Ef 2 < 1, where E is the energy of the particle and 
e <C 1 is the ratio between the straight line and the circle 
radius Q. The billiard dynamics is well described || in 
terms of the map (Q) via the substitutions k = 2eVE, 
T = y/2/E. It is then easy to verify that the quantum- 
integrable regime found in Ref. |l6| E^ 1 < e < i?~ 1//4 , 
coincides with the regime dominated by classical Cantori 
T < k < 1/VT. 

This may be a first indication that not only the classi- 
cal, but also the quantum dynamics of the Stadium, can 
be described in terms of maps: this will be the subject 
of a future work JH| . 

The author is thankful to G.Casati, I.Guarneri and 
D.L.Shepelyansky for useful discussions. 



The dynamical localization mechanism is not con- 
nected with this power law decay. Indeed the same alge- 
braic decay can be found for any k, in the region k < k cr 
as well for k > k q i (at least in the tails of the distribu- 
tion S). In more details, on semiclassically approaching 
the border k q i , the quantum distribution shows big peaks 
of probability for high momentum values which indicate 
that new regions of the classical phase space are now 
quantally accessible. It is exactly the presence of such 
peaks which causes the large increase of l a . The pres- 
ence of bumps of probability far from the initial state 
no = is shown in Fig. 4b. 

In conclusion, a discontinuous map which is a simple 
generalization of the Chirikov standard map has been 
studied. Differently from the latter, the dynamics is 
slowly diffusive even when the motion described by the 
CSM is prevalently regular. In this region the quantum 
analysis reveals quite unexpected features. Above the 
Shuryak border k > T, two different scaling laws for the 
localization length are found. The first, l a ~ y/k/T, 
marked by the presence of classical Cantori acting as to- 
tal barriers to quantum motion, is a region of quantum 
integrability. The second is a region characterized by 
dynamical localization (l a ~ D) thus indicating the exis- 
tence of this phenomenon even in case of slow diffusion. 
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